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ABSTRACT 

The presence of convective and turbulent motions, and the evolution of mag- 
netic fields give rise to existence of temperature fluctuations in stellar atmo- 
spheres, active galactic nuclei and other cosmic objects. We observe the time 
and surface averaged radiation fluxes from these objects. These fluxes depend on 
both the mean temperature and averaged temperature fluctuations. The usual 
photosphere models do not take into account the temperature fluctuations and 
use only the distribution of the mean temperature into surface layers of stars. We 
investigate how the temperature fluctuations change the spectra in continuum as- 
suming that the degree of fluctuations (the ratio of mean temperature fluctuation 
to the mean temperature) is small. We suggest the procedure of calculation of 
continuum spectra, which takes into account the temperature fluctuations. As a 
first step one uses the usual model of a photosphere without fluctuations. The 
observed spectrum is presented as a part depending on mean temperature and 
the additional part proportional to quadratic value of fluctuation degree. It is 
shown that for some forms of absorption factor the additional part in Wien's 
region of spectrum can be evaluated directly from observed spectrum. This part 
depends on the first and second wavelength derivatives, which can be calculated 
numerically from the observed spectrum. Our estimates show that the tem- 
perature dependence of absorption factors is very important by calculation of 
continuum spectra corrections. As the examples we present the estimates for a 
few stars from Pulkovo spectrophotometric catalog and for the Sun. The influ- 
ence of temperature fluctuations on color indices of observed cosmic objects is 
also investigated. 

Subject headings: Radiative transfer, turbulence, temperature fluctuations; Stars: 
atmospheres, color indices; Active galactic nuclei 
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Introduction 



From investigation of the Sun (see lStix Ill99ll ) we know that the photosphere has chaotic 



changes of the temperature. They are due to convective and turbulent motions, and the 
existence of inhomogeneous magnetic fields, which are the specific examples of stochastic 
media. Recall that stochastic behavior of gas motions arises due to many types of instabilities 
occurring in these media. Most important is the convective instability which results in the 
star cells of rising and descending motions of hot and cold gas, correspondingly. For the Sun 
the ratio of mean temperature fluctuations T' to the mean temperature (T) = Tq is equal 
to ?7 = a/ (T'^) /To ~ 0.03. It appears this estimation is minimal because the solar spots are 
not taken into account. No doubt, the temperature fluctuations exist in many stars, not 
only in the Sun. The Sun is the only star which can be observed detailed in every point 
of the surface. Further we restrict ourselves by consideration of point-like cosmic objects. 
Apparently the most high fluctuations exist in stars with the profound convective zone and 
in active galactic nuclei. 

In stochastic media all the values - the temperature T = {T)+T', the radiation intensity 
/ = (/) + the absorbtion coefficient a\ = {ax) + a'^^ and so on, are stochastic, i.e. 
are characterized by its mean values, for example, (T) = Tq and by fluctuations T'. The 
mean value of fluctuations are equal to zero (for example, (T') = 0). The mean values are 
determined by the average over ensemble of realizations during the time of observation. The 
chaotic events in space and time with characteristic life-time greater than the observation 



time are taken into account in the mean values (see, in more detail, iLevshakov &: Kegel 



19971 ). From the stars and other point-like cosmic objects we observe the mean radiative 
flux (H). We stress that the influence of temperature fluctuations on the spectra is pure 
statistical effect, depending on the specific ensemble of temperature realizations in space and 
time, and on the time of observation. It is known that the averaged spec tra of stars, even 



neigh boring in spectral classes, frequently differ one another (see catalog lAlekseeva et al. 



19961 ). The possible reasons may be different, and the influence of temperature fluctuations. 



without doubt, is one of the most important between them. 

In present paper we investigate the influence of temperature fluctuations on the con- 
tinuum spectra of various cosmic objects. The preliminary consideration of this subject is 



given in the paper ( ISilant'ev &: Alexeeva II2008I ). Their paper deals only with investigation 



of effective thermal sources in stochastic media Sx = {ax(T)BxiT)) as compared with the 
usual sources, ax(To)Bx(To), taken without the influence of temperature fluctuations. 

Here Bx{T) is the intensity of Planck's radiation. It was shown that the source Sx can 
fairly strong differ from ax{TQ)Bx{TQ), especially at high value of parameter hu/kTQ. The 
value Sx can be both higher than ax(TQ)Bx(TQ) and lower one, depending on speciflc form 
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of absorbtion coefficient ax{T). Below we will obtain the formula for corrections to observed 
radiation flux, which arises due to existence of temperature fluctuations. 

Before start a detail consideration let us show in a simple example of averaging over two 
temperature reahzations {T — Tq + T' and T — T^ — T'), how arises this statistical effect. 
The average value of Planck's function (in the Wien hmit) is equal to: 



Here we used the condition T'/Tq <^ 1. 

It is seen that the observing averaged value {B\{T)) is higher than B\{Tq). At T'/Tq = 
0.05 and hu/kTo = 10 this growing consists of 13%. At the Rayleigh limit {hu/kT <^ 1) 
the Planck function Bx ^ T and {Bx) = Bx{Tq). For intermediate cases one has {Bx{T)) > 



The analogous calculations demonstrate that the mean value of absorbtion factor {ax{T)) 
can be both higher than ^^(To) and lower than this value, depending on the specific form of 
ax{T). Joint consideration of Planck function and absorbtion coefficient can give rise both 
to positive and negative additional contribution to the usual spectra, depending on the mean 
temperature Tq. 

Below we investigate the problem how the temperature fiuctuations change the con- 
tinuum spectra of cosmic objects in detail. For a number of stars and wavelengths we 
derived numerically the corresponding corrections to the spectra, depending on the degree 
of temperature fluctuations. For the cases, when the main contribution is due to bound-free 
transitions in negative hydrogen ions, we suggest the new techique to calculate the correc- 
tions to the spectra directly from the shape of observing spectra, without using the results 
of any photosphere models. Note that we consider a degree of temperature fluctuations rj as 
a given parameter. The inverse problem - the estimation of 1] from observed spectra is not 
considered in our paper. We hope investigate this difficult problem in the future. 






Bx{To). 
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2. Basic formulas 



see 



The radiation flux in continuum H^^ is related with the Planck function Bx as follows 



Sobolevlll969r ): 
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The functions E^ir) obey the relations: 



(2) 



nEn+i{T) = e 
dr 



- rEn{T), 
-En-iir). 



(3) 



At high T on e exists Eni r) — )■ exp(— r)/r. The flgures Ei{t),E2{t) and E^i^r), for example, 
are given in (IGray 1 119761 ) . 



In Eq.(l) the mean temperature Tq is considered as known function of the optical depth 
Tx. To flnd the relation between the mean temperature Tq and optical depth tx is the main 
task at construction of photosphere models. Formula (1) does not take into account the 
influence of temperature fluctuations. 



I\ \t\,ij) was derived. This eq uation has the effective source Sx, which was investigated 



In the paper ( jSilant'ev II2005I ) the radiative transfer equation for mean intensity (/a) 



by ( jSilant'ev fc Alexeeva II2008I ). The sign () in radiative transfer equation (4) means the 
averaging over ensemble of various realizations along the wave path, i.e. one takes into 
account the local temperature fluctuations. 

In LTE approximation the transfer equation for I^^\tx, acquires the form: 



(nV)/r 



{ax)I^^^ + Sx{To,v), 



(4) 



SxiTo,v) = {ax{T)Bx{T)). 



(5) 
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As a result of averaging, all the values acquire the dependence on the degree of temperature 
fluctuations: 

-to 

The fluctuations T' have the local character. Thus, in the case of the Sun they describes 
the difference between hot and cold parts of granula. The transfer equation (4) is station- 
ary and differs from the usual transfer equation only by the feature that the temperature 
fluctuations are taken into account in absorbtion coefficient ax{T) = ax{TQ + T'). Usually 
in photosphere models one assumes that all physical values are distributed sperically sym- 
metric. We assume the same for temperature fluctuations, i.e. the degree of fluctuations rj 
is independent of the place in visual semisphere of a star. 



According to (jSilant'ev 1120051 ). more exact expression takes into account the correlation 
length R of temperature fluctuations, i.e. instead of term {ax) one uses the expression 
(ax) ■ (1 + {{a'Y)/{{a)Y ■ tx). Here tx = (ax) ■ R is the mean optical depth of characteristic 
length of a turbulence R. Usually one exists tx « 1, i.e. the additional term is of the 
third degree of smallness and can be omitted. Such case takes place for solar granulations 
with R ^ 1000 km in the region of wavelength, where the basic absorbtion is determined by 
negative hydrogen ion. 

From Eqs. (4) and (5) one follows the expression for observed radiation flux {Hx), which 
takes into consideration the influence of temperature fluctuations: 



{Hx{r^)) = 2'n T rf/i/i/f (0, /.) 

^0 



27r / d{rx)E2{{Tx)) — . 7 

Jo MT)) 

Mean values (q;a(T)) are derived by the gaussian formula (see Silant'ev & Alekseeva 
2008): 

{ax{T)) = ^^ dxexp(-—)ax{To-{l + x)). (8) 



^1] J-oo V 
This expression is valid at t] < 0.2. 

When the fluctuations are absent, rj = 0, formula (7) reverts to usual expression (1) 
for where drx = ax{To)dz and Tq = Tq{tx). For brevity we will use simply Tq. Recall 
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that the radiation flux without consideration of fluctuations we denote as . The value 
{ax{T)Bx{T)) is calculated from Eq.(8) analogously. 

Our next task is to relate the observing flux {Hx{7])) with the i.e. to show how 

the small temperature fluctuations change the continuum spectra. As an initial approxi- 
mation we take the expression and then we will seek the corrections to this spectrum 
proportional to rj"^. In section 4 we will consider this problem in Wien's region of spectrum, 
and for such temperatures and wavelengths where the main absorbtion is due to bound-free 
transitions of outer electron in negative hydrogen ion. In this case the flux {H) can be found 
by numerical differentiation of observed spectra over wavelength A. 

The values of fluctuation degree rj are determined by physical conditions in every specific 
object. We assume r/ as a known parameter and derive the formula for estimates of spectra 
corrections due to temperature fluctuations for a number of objects, where absorbtion coef- 
ficient deals with bound-free transitions in the negative hydrogen ions. Of course, it would 
be very interesting and important to derive the technique how to estimate the parameter 
77 from the analysis of observed spectra. It appears this complex task may be resolved by 
statistical comparison of objects with very close physical conditions. In this paper we do not 
consider this difficult problem. 



Assuming that the temperature fluctuations are low (77^ <^ 1), let us use in Eq.(8) the 
expansion series: 



and take into mind that (x) — and (x^) — rf'. Dimensionless parameters 6a (Tq) and a\{TQ) 
have the forms: 



3. Relation of flux {Hx{r])) with Hi 



ax{To ■ (1 + x)) = ax{To){l + bx{To)x + ax{To)x^ + ...) 



(9) 



bx{To) 



To dax{To) 



ax{To) 




a A (To) 



(10) 



Thus, the mean absorbtion factor acquires the form: 



{ax)^ax{To){l + ax{To) rf). 



(11) 
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The optical depth {t\) in expression (7) is calculated according to formula: 



{tx)= dz{ax{z)) =tx + 7]^ I dr'^axiU), (12) 
Jo Jo 

where we use the same notion drx = ax{To{z))dz, as in calculation of the flux H^^ (see Eq. 

(I))- 

Calculating by analogy the value {ax{T)Bx{T)) , we obtain from Eq. (7) the following 
expression for observed flux {Hx{ri)): 



2nr]^ / drxE^irx) 



dBxjTo 



+ 



rp2 



d^Bx{T,) 



(13) 



Ax = 2Tirf drxBxiTo 



The value Tq 



E2{rx)ax{To)-Ei{rx) I dr'^axiU). (14) 

To(rj^). The negative term with axiTo) in Eq.(14) arises from the expansion 
of the exponent exp(— (ta)) (see formula (12)) in power series with parameter rf'. 

Using the known photosphere models, i.e. the dependence between To and optical depth 
Ta, one can calculate the right part of Eq.(13) and obtain the flux {Hx{ri)), corresponding 
to level of temperature fluctuation rj. The second and third terms in Eq. (13) characterize 
correction to model spectrum H^'^ due to temperature fluctuations. In principle, the condi- 
tion of the best coincidence of expression {Hx{ri)) with observed flux allows us to estimate 
parameter r] in the atmosphere of a cosmic object. The influence of fluctuations is most 
strong in the region of short waves. It appears we have to became the best coincidence just 
in this wave region. Taking different values for parameter rj, one can obtain more satisfactory 
coincidence of flux {Hx{r])) with the observed spectrum than that for model case H 



(0) 



The authors are not specialists in the model construction activity. Therefore we restrict 



ourselves by the estimate of difference {{Hx{r])) — H^^^)/H^' only for a number of Sun-like 
stars and the Sun for wavelengths where main contribution to opacity gives the bound- 
free electron transitions in negativ e hydrogen ions if ~. These estimates follow directly 
from observed spectra (see catalog lAlekseeva et al. Ill996l ). For these wavelengths we can 
substitute the temperature derivatives by numerical derivatives over wavelengths in observed 
spectra (the results are presented in Tables 2 and 3). 

First we yield this substitution for Plank's function Bx{Tq). Then such substitution will 
be made for some types of absorbtion coefficients. Note that not all types of a(To) allow us 



r(0) 
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such substitution. Finally, we will consider the spectrum in Wien's region, simultaneously 
using the specific form of bound-free transitions in negative hydrogen ion H~ for those cases, 
where this absorbtion plays the dominant role. 



3.1. Substitution temperature derivatives by wavelength derivatives 

The Planck function has the following form: 



2hc^ 1 
B^{T) = —-f{\T), fiXT) 



^0 _ hv 14388 



(15) 



AT kT XT 

Here and what follows the wavelengths always are taken in microns. The function /(AT), 
depending on the product A and T, obeys the relations: 

These relations allow us to substitute the derivatives over temperature by the derivatives 
over wavelength. The corresponding formulas are the follows: 



dT TV dX 



i,(20i?, + 10A^ + A^^). (17) 



dT^ T2 V dX dX^ 

Absorbtion coefficients ax in c ontinuum sp ectra often have the vast regions of smooth 



dependence on the wavelength (see iGray Ill976l ). In these regions they can be approximated 



by power law a\ ~ a(T)A". Then the values 6a(To) and axiTo) (see Eq.(lO)) are independent 
of A: 6a (To) = &(To) and aA(To) = a(To). In some cases the absorbtion coefficient has the 
form of product q;o(A) and ai(T), i.e. axiT) = ao{X)ai(T). Very important example of such 
dependence is the case of bound-free transitions in negative hydrogen ion H~ (see in more 
detail below). In such cases the values ax{T) and bx(T) also are independent of wavelength. 

Using the relations (17), Eq.(13) can be written in the form: 

{Hxiv)) = Hfh 
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21X7] 



d\ 2 dX^ 



2nri'' I dTxE2iTx)b{To) (^55a(To) + A^ ] + A^. 



;i8) 



The second term in the right part of this expression describes the contribution of tem- 
perature fluctuation to the spectrum due to influence of Planck's function, i.e. this is the 
contribution of value (B \). The estimate of this contribution from model spectra, for ex- 
ample, from models of (IKurucz et al. Ill974j ). is of interest for knowledge of influence of 
temperature fluctuations on continuum spectra. Beforehand, note that the contribution of 
{B\) is smaller than that of term with 6(To). 

Because tx and E2{tx) depend on wavelength A, the A - differentiation cannot be taken 
out the integral in Eq.(18). But for absorbtion coefficients, mentioned above, exist the 
relations: 



-g^ = P(A)aA, -g^ = p{X)Tx, 



/3(A) 



1 dao{X) 
ao{X) dX 



which lead to the following formulas: 



(9R /^tjC) 
27rA / drxE^irx)^ = X^^ - X(3{X)Tx, 



dX 



dX 



(19) 



Tx = 2Hf - 27r/r(0, 1) = 4n I rf/i/i[/r (0, /i) - /r(0, 1)] 



r(0) 



r(0). 



r(0), 



(20) 



27rX^ / drxE^irx) 



d'Bx 
dX^ 



(0) 



aA^ 



_ dm 



dX 



Tx 



f)T 

2A^/3(A)— ^ + A2/32(A)Ma. 
aA 



(21) 



When deriving these formulas, we used Eqs.(l) and (3). The values /3(A) and d^/dX for 
negative hydrogen ion are present ed in Table 1. They were calculated numerically from 
figure for ax, given in ( iGray 1119761 ). 



Function Tx depends on the flux of outgoing radiation H^^ and the intensity of radia- 
tion /f^(0,l) escaping perpendicular to the surface. Note that due to the law of radiation 
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darkening {lf\o,fi) < lf\o,l)), the function Tx is negative. This function can be calcu- 
lated from a model of photosphere. Here we will use the simplest model, assuming that the 
function Bx{tx) ii icreases linearl y with r^. This model is presented in a number of books 
(see, for example, ISobolev 1119691 ). 



Below we restrict ourselves by the Wien limit of a spectrum, when the parameter 
Qq/XT » 1. In this case simple model, mentioned above, gives rise to the relations: 



-H 



(0) 
A ' 



Mx = 2Tx + 2n / dTxBx{To){Tx - 1) exp (-r^) ~ 



(0) 
A • 



(221 



According to formulas (19) and (21), contribution of {Bx{tx)) into {Hx{ri)) acquires the 



form: 



{Hxiv)) = 



(0) 



(0) 



2 



A^/3(A) ^ - (^5A/3(A) + 3^'^ ) + -A^/3^(A)M, 



(23) 



This formula depends on the specific form of absorbtion coefficient. Because of the difference 
{Bx{Tq, 7])) — Bx{Tq) is positive and slowly decreases with grow of wavelength, then the term 
{Bx) at rj"^ gives the positive contribution to observed spectrum {Hx{ri)). Yet the contribution 
of term bx(To), in principle, can exceed the contribution from {Bx{To)), and, in particular, 
gives rise to negative correction to (i^A)- 



3.2. Estimate of expressions with ax(To) 

The term Ax in formula (14) describes the change of spectrum as a result of difference 
between (ax) and ax{To) (see formula (11)). Negative part of Ax presents the decrease of 
flux due to additional term with 77^ in exp(— (ta)) (see (12)). The positive part describes the 
increasing of flux due to effective source (axBx). The difference of these terms is lower than 
every separate term. 

It is of interest to estimate the contribution of term Ax to expressions (13) and (18). 
Let us replace the terms axiTo) and ax{TQ) by effective mean value ax, taken at the level 
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Tx — 2/3, where is located the main source of radiation. As a result, one can obtain the 
expression: 



POO 

Ax - 2TrT]^ax / drx {2E2iTx) - e'^') Bx{T^) 
Jo 



2717] ax / dT\TxE2{Tx 



dBxin 



rfaxTx ^ -r]'^ax6xH 



2— 



drx 
r(0) 



(24) 



This expression is less than the term lOH^^^ in Eq. (23) if the value ax is less or of the 
order of unity. In these cases the term Ax can be neglected. The coefficient 6x, according 
to approximate Sobolev's (1969) formula, is estimated as the ratio of the mean absorbtion 
factor over all spectrum to the local coefficient ax- I n the case of absorbt ion by negative 
hydrogen ions H~ the Chandrasekhar's estimation (see lChandrasekhar 1119501 ) gives the value 
6x — 0.26. This gives the estimation Ax — 
formula (23). 



1 2)Hf'\ which is less than the term lOHf'^ in 



4. Estimates in Wien's limit 

The case go/{\TQ) = 14388/(ATo) ^ 1 corresponds to Wien's limit of thermal emission. 
Recall that at the maximum of Planck's emission XmaxT = 2900. This gives Qq/ (XmaxT) — 5, 
i.e. the Wien approximation Bx{T) = {2h(? /\^) exp [—go/XT) exists at least up to maximum 
of spectrum. The relative deviation from Planck's formula at maximum of emission consists 
of 0.7%. At go/ XT = 3.6 this deviation is ~ 2.8 %, that corresponds to the value X{fim)T ~ 
4000. 



This means that the Wien limit occurs practically up to A = 0.5 /zm at T = 8000 K, 
and up to A ~ 1 /im at T ~ 8000 K. Note that in this temperature interval the absorbtion 
coefficient in continuum, ba sically, is due to bound - free transitions of outer electron in neg- 



ative hydrogen ion H (see IChandrasekhar I Il950l ) . Below we will give the formulas of this 



case, in detail. In Wien's limit the Einshtein-Milne correction due to stimulated emission 



(the factor (1 — exp {—go/ XT)) is low and the formulas for absorbtion coefficient (see I Gray 



19761 ) become simpler. 

For the absorbtion factors of the type 



ax = ao(A)T^e-^i/^ 



(25) 
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dimensionless coefficients b{T) and a{T), according to Eq. (10), acquire the form: 



a(T) = 2(1^ + (izik + 



T 



2T2- 



(26) 



Formula (25) is valid for the description bound-free transitions in H , where 7 = —1.5 
and gi — —8700 (if the wavelengths are in microns). 

The coefficients b{T) and a(T) obey power-law dependence on the temperature. There- 
fore in Wien's approximation Eq.(18) (without the term Ax) can be simplified if we transform 
the terms of type Bx{Tq)/Tq to expressions with the differentiation over wavelength. For 
this aim we can use the following equalities (introducing the notion fi{XT) — exp{—gi/XT)): 



/i(AT) 9/1 /i(AT) 



2 a/i d'h 

X dX dX^ 



(27) 



T dX' T2 v^i. 

Taking into account, that in Wien's approximation Bx ~ A~^ exp(— g^o/AT), we derive the 
relations: 



.dBx 
dT 

Bx 



A2 
9l 



30Bx + 12A 



dBx , ,,d'Bx 



dX 
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Using these relations, and Eqs. (19) and (21), we derive Eq. (18) (without the term 
proportional to a{T)) in Wien' approximation: 

{Hxiv)) = iff + 



Co{X)Hf + Do{X)Tx + Go{X)Mx + Ci(A)A 



Z..(A)A^ + C.(A)A^^ 



dH 



(0) 



dX 



(29) 



Co(A) = 10 + 57 + 30e, Ci(A) =5 + 7 + 12^, 
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C2(A) = (0.5 + O, 



DoiX) = -A/3(A)Ci 



GoiX) = yi3\X)C2, A(A) = -2A/3(A)C2 



(30) 



Here we use the dimensionless values ^ = QiX/qq with qq = 14388. For the case a\{H^) 
one has 7 = —3/2, gi = —8700, gi/go = —0.6 and ^ = — 0.6A. As it was mentioned earher, 
we can negle ct the term A\ if the absorbtion in negative hydrogen ions are considered. 
According to ( IChandrasekhar Ill950l ; ISobolev Ill969l ) this case is valid for the solar type stars 
with Te < 7 ■ 10^ K. 

In Table 2 we present the corrections to {Hx{r])) for a number o f solar type stars (/3 Hy i 
(G2IV), /i Vel (G5III), l Per (GOV)) and the Sun, using the catalog JAlekseeva et al. Ill996h . 
The corrections were calculated for A = 0.55/im and 0.6 fim. For these wavelengths one can 
neglect the contribution of free- free transitions in H~ . Note that the columns with positive 

and the columns with negative numbers present the values 



numbers in Table 2 refer to 



$A- The notions of these values are given below. 

Formula (29), with taking into account Eq.(22), acquires the form: 



^2 H-(O) 



dX 



dX^ 



(31) 



where we introduce the notion (3' = 8(3/ dX, and the factors Aq and Ai have the following 
form: 



Ao = (5 + A/3)7 + e(30 + 12A/3 + X\p^ + /3')), 



= 7 + e(12 + 2A/3). 



(32) 



The term coincides with that in square brackets in Eq. (23), if the relations (22) 

are used. This term has the form: 



^(0)^(0) 



A^ 



[10 + 5A/3 + y(/3^ + /3')]i^A 



(0) 
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+(5 + A/3)A ^ 



A2 S^i/f 



2 

Recall that this term presents contribution of (-Ba)- 



9A2 



(33) 



4.1. Results of calculations and the discussions 

Usually the star fluxes are presented in logarithmic scale Igi^A or in star magnitudes 
mx = —2.5lg H\. Our formulas have the derivatives dH\/d\ and d^Hx/ . Differentiation 
of the relation In Hx = In 10 Ig Hx gives rise to expressions: 

-9^ = H.\nlO^^ = -Hx^^, (34) 



9A2 



mio^+rinio^^^^^ 



9A2 



dX 



In 10 drri). 



2.5 9A 



In 10 d'^mx 
T5" 9A2 



(35) 



These formulas show that the derivatives are proportional to Hx- It means that the 
expression (31) may be presented in the form: 



H 



(0) 



H 



(0) 



(36) 



The value AHx can be presented as a corresponding change AT (see ISobolev Ill969l ). As a 
result, we derive the estimate: 



AT ^ (2.8 - 4) (^^^ r/2^,To, (37) 

Instead of Tq in such estimates one may be taken the effective temperature Tg in the regions, 
where BxiTe) satisfactory describes the spectrum. Recall that in Wien's region even low 
temperature changes AT give rise to fairly high changes of fluxes. 

Let us present the explicit formula for function $a, if the absorbtion coefficient has form 
(25). This formula follows directly from Eqs. (31) - (35): 
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10 + 5A/3 + y(/32 + /3')+Ao+ 



(5 + A/3 + A In 10 



d}gHx 
dX 



+ 0.5) A^ 



mio^ + Anio^^^^^ 



dX 



(38) 



First of all, we obtain the part of expression ^x, which is due to contribution of mean 
Planck's function (Bx). This expression can be obtained from Eq. (38), where terms with 7 
and ^ have to be omitted. Let us denote this expression as ^^x ^ ■ To calculate this term one 
can use the model spectra, or use the observed spectra {Hx{ri)), assuming that at low values 
T] exists approximate equality {Hx{ri)) ~ Hx^ ■ 

As it was mentioned earlier, the spectra depend strongly on the contribution of the 
coefficient 6a(7o) (see Eq. (10)), which describe the change of effective source (axBx) due 
to fluctuations of first temperature derivatives in absorbtion factor ax{T) and in Planck's 
function Bx{T). As a role, the first derivatives characterize the changes most sensitively than 
the second ones. The contribution of second order derivatives give rise to less contribution 
to changes in spectra. For this reason, the calculation of ^x^ is useful for comparison with 
basic terms proportional to mean value of the product of first order derivatives. 

It is of interest, that calculations, presented in Table 2, show that the terms with /3 
and /?' give comparatively low contribution to $a , less than (2-^5)% compared with other 
terms. Note, that this low contribution is not the consequence of smallness of j3 and /?' (see 
Table 1 ) , and basically are due to powers of small wavelengths (A ~ 0.5 /xm) at coefficients 
(3 and /?'. Therefore for estimations of {Bx) (term $^°^) can be used the simple formula: 



^(0)^(0) 



lOH^^^ + 5A 



dH 



(0) 



A2 d^Hf 



dX + 2 9A2 • ^^^^ 
The term with bx can be simplified accordingly. The sense of these simplifications consists 
of conclusion, that to obtain the estimates in Eq. (18), one can take out of the integrals the 
differentiations over wavelength. Apparently, this can be used in other cases, if the absorbtion 
coefficients change slowly in the range of considering spectrum, and if the spectrum itself 
has fairly smooth form. In this case one can introduce the large-scale differentiation over 
wavelengths for spectrum itself, and use the averaged absorbtion factor. 

In Table 3 we presented the valu es of functions $1°^ and ^x ^oi stars a PsA and /3 Ari, 



calculated from homogeneous catalog (lAlekseeva et al. Ill996l ). The tables of star magnitudes 
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mx in this catalog are presented with the interval AA = 0.0025 fim. The first star is of the 
class A3V with ~ 8300 K, and the second one is of the class A5V with ~ 8500 K. 
Spectra of these stars does not determined by absorbtion in negative ion H~ . To receive 
the esimates of tem perature fiuc tuations infiuences we used the absorbtion coefficients in 



continuum given in (lAUen 1 119731 ) . We used the parabolic interpolation of data from this 
book for intervals A = 0.4,0.5,0.8 /im and for 6 = 5040/T = 0.4,0.5,0.6 at the value of 
electron pressure IgPg = 1- As a result, we derived the formula: 

a,(e) = 10^(^'®), 
C(A, 0) = —[\\20e^ - 210 + 5.3) + 

o 

A(-2402 + 25.50 - 6.33) + (5.80^ - 6.180 + 1.468)]. (40) 

Taking into account that A - dependence in this expression is fairly smooth, we used for 
estimations the mentioned above (see Eq. (39) and so on) simple expressions. 

The negative corrections $a denote that the temperature (see the notions (36) and 
(37)) is less compared with the temperature obtained from photosphere models, which do 
not take into account the temperature fiuctuations. Especially important conclusion from 
our estimates is that most considerable corrections to {Hx{r])) are due to terms with bx(T) 
(see Eq. (10)). These terms give rise to negative corrections whereas the contribution of 
{Bx{T)) gives positive corrections. The contribution of terms with 6a is determined by the 
average of product of first temperature derivative from absorbtion factor and from the Planck 
function, i.e. most important linear expansions over T' are averaged. 

It is of interest to note that observed Sun's spect rum in ultraviolet region is lower tha n 



that determined by accepted effective temperature (jStix I Il99ll : IChance &: Kurucz Il2010l ). 
Usually this effect is related with influence of absorbtion of metals. However, the possible 
contribution of temperature fluctuations is also not eliminated. Though this mechanism is 
not studied in detail, but preliminary we can note that the calculations of source function 
of new transfer equation (4) in ultraviolet region demonstrate the strong difference from the 
usual Planck's function. 

The estimates (37) for the Sun at A = 0.55 /im and Te ~ 5770 K give AT ~ -(2.8^4)7/^- 
9.8re ~ -(2.8^4) -12470772, i.e. AT ~ -(31^45) K for = 0.03, and AT ~ -(87^125) K for 
T] = 0.05. In percents these estimates are (0.5-f-0.8)% and (1.5^2.1)% of Tg, correspondingly. 

These estimates are given for A ~ 0.55 fim and practically do not depend on low changes 
in taking Tq ~ Tg. Note, that to estimate the change in effective temperature of a star, one 
has calculate the change of all the photosphere spectrum in the range of all wavelengths. 
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5. Influence of temperature fluctuations on color indices of radiating objects 

The color index of observing objects is characterized by the difference of radiation 
fluxes between standard wavelength intervals. Usually the observed ffuxes F\ correspond to 
standard wavelength intervals near the central wavelengths (Aj/ = 0.36/im, = 0.44/i,m, 
Ay = 0.55yU.m, = 0.7 fim). For example, the color index U — B is determined according to 
formula: 



p(0) 

[/(o) _ ^(0) _ _2_5 Ig + Const. (41) 

The value Const is taken from the reasons of convenience. Frequently one takes Const — 
for stars of AOV-type. 

Usually one considers that fluxes F^'' = Bx{To), where temperature Tq on an average 
describes the observed part of spectrum. This relation formally follows from usual transfer 
equation at supposition that the source of radiation Bx{Tq) is independent of optical depth 
in radiating region. The radiating objects can have the fluctuating temperature. In this 
case the source is described by the formula {ax{T)Bx{T)) /{Bx{T)) = Sx{To, rj) /{Bx{T)) (see 
Eqs. (4) and (5)). 

As we know, in this case the flux {Hx{r])) = (-Fa) depends on the degree of temperature 
fluctuations rj. Assuming that this effective source is independent of optical depth, we obtain 
the expression (Fx) — Sx{To,r])/{Bx{T)). This means that the color indices (U) — (B) and 
so on depend on degree of ffuctuations f]. As usually, we restrict ourselves to low ffuctuations 
<^ 1 and instead of Eq. (41) obtain the relation: 

(U) - (B) = C/(°) - + r]^A{U - B), (42) 
where A{U — B) is expressed by the formula: 



A{U -B) = -2.5{Kx^ - K, 



K ^0 d'BxjTo) To dBx{To) 

' 2B,(To) dTi ^ Bx{To) dTo ' °^ 

= xf + X«6,(To). (43) 

Analogous formulas characterize the color indices B — V and V — R. The numerical values for 
coefficients K^x\'^o) and K^x\'^o) presented in Table 4 for wavelengths A = 0.36; 0.44; 0.55 



- 18 - 



and 0.7 fim. Note, that these factors are positive. At growing of wavelengths these coefficients 
diminish. 

It is seen from Eq.(43) that the color indices strongly depend on specific form of ab- 
sorbtion coefficients. This is seen clearly from the example, when we use the absorbtion 
factor (25) for bound - free transitions in negative hydrogen ions. In this case b\{To) has 
the form (26) with 7 = —1.5 and gi = —8700. The term K^^ describes the influence of 
temperature fluctuations without contribution of absorbtion factor. This term is positive, 
as it was mentioned earlier. The expressions A{U — B),A{B — V), and A{V — R) with 
allowance for only K^^^ - term are negative for all mean temperatures Tq. The consideration 
of absorbtion coefficients gives rise to the result that the total coefficient K^^ + A'^^^ &a(7o) 
strongly differs from the term K^*^ (see Table 5). In Table 5 we also present the values of 
A{U — B) and so on, which take into account both coefficients - K^^ and 

At comparatively low temperatures the correction to color indices is negative. Then, 
with the grow of the temperature, this correction acquires positive sign. Most high the 
influence of temperature fluctuation is in color index U — B, i.e. for small wavelengths. It 
seems, that in active galactic nuclei the temperature fluctuations may be high. In these 
cases the consideration of fluctuations can change the color indices considerably. Besides, 
at introduction of color indices would be useful, though roughly, take into account the 
temperature gradient nea r the surface of emitting object, for example, using the mentioned 



approximate formula in ( ISobolev I Il969l ). It appears the inverse problem of estimation of 
degree of temperature fluctuations rj is simpler from consideration of color indices of cosmic 
objects with like physical conditions. 



6. Conclusion 

Let us present short discussion of obtained results. The existence of stochastic processes 
in photospheres of stars, active galactic nuclei etc., which arises due to convective and tur- 
bulent motions, generation and evolution of magnetic fields, give rise to stochastic behavior 
of temperature, intensity of radiation and absorbtion coefficients. These stochastic processes 
determine the mean radiation flux (Hx). This flux differs from the flux which takes 

into account only dependence of all values from mean temperature Tq and is independent 
of temperature fluctuations T'. The numerous models of photospheres determine the mean 
temperature as a function of distance from the photosphere's surface. 

Assuming the degree of temperature fluctuations rj = (T'"^) /Tq as a low value, we ob- 
tain the relation between the observed flux {H\) and model flux H^^'' ■ difference between 
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these fluxes is proportional to yf and can assume both positive and negative values. This 
is determined by the specific form of absorbtion coefficient. In Wien's hmit this difference 
can be estimated from the observed spectrum, because it depends of the first and the second 
derivatives over wavelength, which can be calculated numerically directly from the observed 
spectrum and is independent of specific choice of photosphere model. 

We calculated this difference of fluxes for the Sun, and solar type stars (/3 Hyi, n Vela and 
i Per) for A = O.SS/xm and A = O.G /zm, using the homogeneous Pulkovo spectrophotometric 



catalog (see lAlekseeva et al. 1119961 ). Also we obtain the corresponding estimates for stars a 
PsA and /3 Ari. 

These estimates demonstrate that the most high correction is due to mean value of 
product of first temperature derivatives of the absorbtion coefficient and the Planck function. 
The estimations show that the difference of spectra (-^a) and is fairly high and can be 
used for correction of model results, especially in ultraviolet region of spectra. We also 
described the influence of temperature fluctuations on the color indices. It appears these 
corrections can be high, especially for active galactic nuclei. 
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Table 1: Functions /3(A) and dp{X)/dX for absorbtion factor ax{H ). 



A, urn 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 
^ 3.5 2.4 1.7 1.28 1 (176 0.58 0.41 0.3 0.22 0.15 

dp/dX -28 -18 -11 -6.8 -5.3 -4.25 -3.5 -2.7 -2.0 -1.5 -1.25 



Table 2: Functions $1 and ^\ for the solar type stars. 







the Sun 




/3Hyi 




// Vel 




L Per 




A = 


0.55//m 


4.3 


-9.8 


3.6 


-23.6 


7.9 


-12.9 


5.9 


-18.0 


A = 




2.7 


-10.2 


3.3 


-26.4 


3.9 


-14.4 


5.6 


-20.1 



Table 3: Functions $3, and for spectra of a PsA and ^ Ari. 



A, /im 


0.45 


0.47 


0.50 


0.55 


0.60 


a PsA, 


1.7 


1.4 


1.0 


0.4 


0.0 


a PsA, $A 


-10.2 


-11.3 


-12.2 


-13.5 


-13.4 


/3 Ari, ^J'^ 


3.1 


2.8 


1.9 


1.4 


0.9 


/3 Ari, *A 


-8.6 


-9.2 


-10.8 


-11.0 


-11.1 



Table 4: The values of K^^ {^0) and K{'>{To). 



To, K 


4-103 


5-10-3 


6-103 


7-103 


8-103 


9-103 


10-103 


12-103 


14-103 


16-103 


18-103 


< 


39.93 


23.98 


15.60 


10.73 


7.71 


5.72 


4.37 


2.73 


1.82 


1.29 


0.95 




9.99 


8.00 


6.67 


5.73 


5.03 


4.49 


4.07 


3.45 


3.03 


2.72 


2.49 


< 


25.27 


14.93 


9.57 


6.51 


4.63 


3.42 


2.60 


1.61 


1.08 


0.77 


0.57 




8.18 


6.55 


5.47 


4.72 


4.16 


3.73 


3.40 


2.92 


2.59 


2.35 


2.17 




14.93 


8.65 


5.46 


3.67 


2.60 


1.91 


1.45 


0.90 


0.61 


0.43 


0.32 




6.55 


5.26 


4.42 


3.83 


3.40 


3.07 


2.82 


2.46 


2.21 


2.03 


1.90 


4'^ 


8.27 


4.70 


2.93 


1.96 


1.38 


1.02 


0.78 


0.49 


0.33 


0.24 


0.18 


5.17 


4.18 


3.54 


3.10 


2.78 


2.54 


2.36 


2.09 


1.91 


1.78 


1.68 
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Table 5: Various functions for bound- free transitions in negative hydrogen ions. 







4-10=^ 


5-10=^ 


6-10=* 


7-10=^ 


8-10=* 


9-10=^ 


10-10=^ 


12-10=^ 






3.21 


-1.92 


-4.08 


-4.98 


-5.31 


-5.36 


-5.28 


-4.96 


Kx, 




-4.79 


-6.29 


-6.57 


-6.43 


-6.12 


-5.79 


-5.46 


-4.87 






-9.14 


-8.40 


-7.57 


-6.83 


-6.20 


-5.67 


-5.24 


-4.56 


Kxa 




-10.73 


-8.84 


-7.51 


-6.54 


-5.81 


-5.25 


-4.81 


-4.16 


A{U- 


B) 


-20.00 


-10.94 


-6.24 


-3.62 


-2.04 


-1.06 


-0.45 


0.21 


A{B- 


V) 


-10.88 


-5.26 


-2.49 


-1.00 


-0.18 


0.29 


0.55 


0.77 


A{V- 


R) 


-3.97 


-1.12 


0.14 


0.71 


0.96 


1.05 


1.07 


1.01 



